The origin of electron-hole asymmetry in graphite 
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The electron hole asymmetry has been measured in natural graphite using magneto-optical ab- 
sorption measurements. A splitting is observed for the transitions at both the Tf-point and the 
//-point of the Brillouin zone of graphite where the effect of trigonal warping vanishes. This result 
is fully consistent with the SWM Hamiltonian providing the free electron kinetic energy terms are 
taken into account. An identical electron-hole asymmetry should be present in graphene. 



The band structure of graphite has been calculated 
by Slonczewski and Weiss (SW) in the late fifties pQ. 
Based upon detailed group theoretical considerations the 
SW Hamiltonian, with its seven tight binding parame- 
ters 70, ..,75, A, can be exactly diagonalized to give the 
band structure. Due to the inter layer coupling the in- 
plane dispersion depends on the momentum k z parallel 
to the c-axis. McClure derived the magnetic Hamiltonian 
for the case when the magnetic field is applied parallel 
to the c-axis [2]. The so called Slonczewski, Weiss and 
McClure (SWM) Hamiltonian has infinite order since the 
trigonal warping term 73 couples Landau levels with or- 
bital quantum number n to Landau levels with quantum 
number n + 3. This coupling breaks the dipole selection 
rule and gives rise to a large number of harmonics in 
the cyclotron resonance. Nakao showed that the infinite 
Hamiltonian can be successfully truncated to a reason- 
able size and numerically diagonalized to find the eigen 
values [3]. At the -H-point the effect of 73 vanishes and 
the SWM Hamiltonian can be analytically solved to give 
a Landau level energy spectrum which depends only on 
7o- This is the origin of a widespread misconception in 
the literature, including our own work, that there is no 
electron-hole asymmetry at the //-point. 

The electronic properties of graphite are well docu- 
mented in the literature [4l418j . In particular, magneto- 
optical techniques have been extensively used to probe 
the Landau level energy spectrum at the H and /C-points 
where there is a joint maximum in the optical density of 
states [TM2~4"] . This data was for the most part analyzed 
using the effective bi-layer model [25] for graphite with 
only two parameters, 70 and an effective inter layer cou- 
pling 27! . The splitting of the /C-point transitions in the 
magneto-reflectance data was analyzed within the effec- 
tive bi-layer model by including electron-hole asymmetry 
due to the non vertical coupling term 74 phenomenologi- 
cally |22j . In our previous work [23j the observed splitting 
of the //-point transitions was not assigned to electron- 
hole asymmetry as there is no trigonal warping at the 



//-point, so the effect of 74 vanishes. 

In this letter we show that electron-hole asymmetry 
exists for all values of k z and is an inherent part of the 
SWM Hamiltonian through the often neglected free elec- 
tron kinetic energy terms. The asymmetry should lead to 
an observable splitting of both the H and ZsT-point optical 
transitions. Extending our previous magneto-optical on 
natural graphite to lower energies, lower temperatures, 
and higher magnetic fields we show that a splitting of 
both the H and /C-point transitions due to the electron- 
hole asymmetry is observed. The size of the splitting 
at the //-point is in good agreement with the predicted 
electron-hole asymmetry. The splitting of the /"T-point 
transitions is also found to be dominated by the free elec- 
tron terms with 74 and 75 playing only a secondary role. 

Nakao [3] derived an explicit form for Landau level 
energy spectrum at the //-point. Unfortunately, when 
writing the expression Nakao neglected for simplicity the 
small free electron kinetic energy terms H 2 k 2 /2m, where 
k is the in plane wave vector and m is the free electron 
mass. The free electron terms are quantized in a mag- 
netic field and their values are significant for all magnetic 
fields. The SWM Hamiltonian can easily be diagonalized 
at the //-point and the correct expression for the Landau 
level spectrum, including the free electron terms is, 

_ A ± ^(A + h 2 s/2m) 2 + 3ns 7 ga§ nh 2 s 



A ± ^(A - fi 2 s/2m) 2 + 3(ti + l)s7o a o (n+l)H 2 s 
E ^ 2 ~ 2 + 2m 

where n = 0,1,2,... is the orbital quantum number, 
s = 2eB/h and ao = 0.246 nm. The Zeeman term 
has been omitted since it simply shifts the energies by 
±gnsB/2 and can easily be added if required. At the 
//-point the electron hole asymmetry is provided by 
the free electron term nh 2 s/2m. Thus the dipole al- 
lowed transitions, E$_ -> E^ 1 and E 7 ^ 1 -> E% + will 



2 



be split by SE = h 2 s/2m ~ 0.23 meV/T. Note that 
h 2 s/2m <C s7oOq so that to a very good approximation 
E^ — E™ 2 i.e. the Landau ladders remain doubly de- 
generate at the //-point. Note that the free electron term 
has the expected phase (0) for massless Dirac fermions. 

In a similar way, the bi-layer expression |25) can be 
modified phenomenologically to include the free electron 
term for massive fermions with a phase of (1/2) at the 
if -point 
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where n = 0,1,2..., A = 2, e = Vfv2e%B is the char- 
acteristic magnetic energy, Vf = \/3eaojo/2h is the 
Fermi velocity. Eq([2]) has not been derived explicitly, 
however, we have verified that the predicted behavior 
is in exact agreement with the SWM calculation with 
73, .., 75, A = 0. The full SWM model has quantum num- 
bers — 1, 0, 1, 2, ... and there are two special Landau levels 
(LL0 and LL-1) whose energy remains close to zero. The 
bilayer model (Eq.Q) correctly predicted the energy of 
LL-1 using n = 0. LL0 is missing but can be reproduced 
accurately between — 150 T using n — if the free elec- 
tron term is replaced by (n + 3/2)fi 2 s/2m- 16(h 2 s/2m) 2 . 

Before presenting the experimental data, the impor- 
tance of the free electron kinetic energy terms is demon- 
strated by numerically diagonalizing the truncated 600 x 
600 SWM matrix for a magnetic field B = 0.3 T using 
the SWM parameters of Nakao [5] to allow a compar- 
ison. The calculated Landau level dispersion along k z 
is shown in Figjlja) including the free electron terms. 
The symbols (circles and triangles) in FigfTJa) are taken 
from the calculations of Nakao at the same magnetic field 
(Fig. 3 of Ref.[3]). The triangles distinguish the triply de- 
generate Landau levels, which have a markedly different 
dispersion along k z and correspond to leg orbits. Clearly 
there is perfect agreement between the two calculations. 
On the other hand, the calculations in Figjljb) which 
neglect the free electron terms are significantly different. 
Notably, the electron cyclotron energy is underestimated, 
while the hole cyclotron energy is overestimated. Thus, 
the free electron terms have to be included in the SWM 
Hamiltonian if the correct energy spectrum is to be ob- 
tained. As our SWM calculations agree perfectly with 
the results of Nakao, we conclude that the free electron 
terms were omitted from Eq(9) of Ref.[3], but included 
in the numerical calculations of Nakao. 

For the magneto-transmission measurements suitably 
thin samples were fabricated by exfoliating natural 
graphite. The measurements were performed in pulsed 
fields < 60 T (~ 400 mS) . A tungsten halogen lamp pro- 




FIG. 1. (Color online) (a) Calculated Landau level dispersion 
(solid lines) along k z using the SWM parameters of Nakao [3j 
and including the free electron terms. For comparison the cal- 
culated values of Nakao (symbols) are shown, (b) Calculated 
Landau level dispersion along k z neglecting the free electron 
terms. 



vides a broad spectrum in the visible and near infra-red 
range and the absorption is measured in the Faraday con- 
figuration with the c-axis of the graphite sample paral- 
lel to magnetic field. A nitrogen cooled InGaAs photo- 
diode array, or an extended InGaAs detector analyzed 
the transmitted light dispersed by a spectrometer. The 
use of two detectors allows us to cover a wide energy 
range 0.6— 1.1 eV. Differential transmission spectra were 
produced by normalizing all the acquired spectra by the 
zero field transmission. Measurements to higher fields 
< 150 T were performed using a semi-destructive tech- 
nique and pulse lengths of ~ 10 fjs and the transmission 
of a polarized CO laser (0.229 eV) measured as a func- 
tion of the magnetic field using a nitrogen cooled HgCdTe 
photodiode coupled with a 200 MHz low noise amplifier 
and an infrared tunable wave plate. 

Representative differential absorption spectra mea- 
sured at T ~ 1.8 K in magnetic fields B — 55 — 59 T are 
shown in Figj2](a) . The spectra contains a large number 
of lines reflecting the large number of K and H point 
transitions which cross in this energy region. Neverthe- 
less, a clear splitting of the ii-point and the if -point 
transitions is observed (arrows). The energy of the ob- 
served transitions are plotted as a function of magnetic 
field in Figj2](b) . Before discussing these results it is use- 
ful to consider the possible transitions at the iT-point. 
Dipole allowed transitions have a change in the orbital 
quantum number of ±1. Due to the doubly degenerate 
Landau level spectrum at the ii-point with E^ 1 = E™ 2 
there are a large number of allowed transitions between 
the valence band (£3- or E2) and the conduction band 
(£34- or i?i). However, the understanding of the prob- 
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FIG. 3. (Color online) (left) Band structure of graphite along 
the H — K — H edge, (right) Schematic of the Landau level 
energies at the //-point showing the electron-hole asymmetry. 
Arrows indicate dipole allowed transitions (An = ±1). Tran- 
sitions are labeled as "effective" E3- — > £3+ transitions with 
"apparent" dipole selection rules An = ±1, 0, ±2. 



FIG. 2. (Color online) (a) Differential magneto-transmission 
spectra of natural graphite measured at magnetic fields in the 
range 55 — 59 T at T ~ 1.8 K. (b) Magnetic field dependence 
of the observed optical transitions in natural graphite. The 
calculated SWM energies of the transitions are shown as lines: 
//-point An = ±1 (thin blue lines), "effective" //-point An = 
±2 (dashed red lines), An = (dotted green lines) and K 
point An = ±1 (thick black lines). 

lem is greatly facilitated by fact that all transitions in- 
volving bands E2 or E\ are degenerate with E^ — > E$ + 
transitions with "apparent" selection rules An = and 
An = ±2. This is shown schematically in Figj3] The 
electron hole asymmetry, also shown schematically here, 
splits both the An = ±1 and the An = ±2 transi- 
tions, while the An = transitions remain unaffected. 
From Eqjl] the splitting of the An = ±2 transitions is 
SE = h 2 s/m i.e. twice the size of the splitting of the 
An = ±1 transitions. 

The energy of the observed H and K-po'mt transitions 
are plotted as a function of magnetic field in Fig[2jb). As 
seen in the raw data, a splitting of the //-point and the 
K -point transitions is observed. The calculated SWM 
transitions energies are indicated by the solid and broken 
lines. 

The energy of the //-point transitions depends only on 
70 = 3.15 eV and the calculated splitting is independent 
of all other SWM parameters and vanishes only if the free 
electron terms are not included in the Hamiltonian. We 



have verified that the predictions of Eq. (nj) are exact. The 
observed splitting of the //-point — > E^t 1 ^ 

transitions (blue solid lines) is beautifully reproduced by 
the calculations. We stress that in either approach there 
are no fitting parameters; the size of the splitting is sim- 
ply given by h 2 s/2m ~ 0.23 meV/T. In addition, the 
observed splitting of the "effective" E3- — > E% + transi- 
tions with "apparent" selection rules An = ±2 (dashed 
red lines) is twice as large in agreement with the predic- 
tions for electron hole asymmetry in Eq([2]). 

The calculated splitting of the /C-point transitions de- 
pends on the SWM parameters used, notably 74 and 75. 
We adjust very slightly 71 = 0.37 eV to fit the observed 
transitions (slope of the magnetic field dependence) and 
use the accepted values for the other SWM parame- 
ters which are summarized in Table ||J The agreement 
turns out to be very good making a further refinement of 
the parameters unnecessary. A comparison of the SWM 
splitting ~ 23 meV at B = 60 T with hs/2m ~ 14 meV 
suggests that 74 and 75 are responsible for approximately 



70 = 3.15 eV 71 = 0.37 eV 72 = -0.0243 eV 
73 = 0.31 eV 74 = 0.07 eV 75 = 0.05 eV 
A = -0.002 

TABLE I. Summary of the SWM parameters used. 
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FIG. 4. (Color online) (a) Magneto-transmission of natural 
graphite showing mainly if-point transitions, (b) Calculated 
SWM transitions together with the measured splitting (sym- 
bols). There is no electron-hole asymmetry for the — > 1 
K-point transition which splits due to the Zeeman term. 



40% of the splitting. The relative importance of the con- 
tribution of the free electron kinetic energy terms to the 
electron-hole asymmetry means that any data analysis 
which neglects them would lead to a significant over es- 
timation of size of 74 or 75 . 

Polarization resolved magneto-transmission, in fields 
up to ±140 T are shown in FigQa). Mainly .PC-point 
transitions are observed in this energy range. The differ- 
ent field directions corresponds to different polarizations 
and the features are shifted in field due to the different 
energy of the n n + 1 and n + 1 — > n transitions. The 
feature around 100 T is the fundamental — >• 1 transi- 
tion, which should not be split (shifted) since the LL0 is 
special and has a free electron term ~ (n + 3/2)h 2 s/2m 
with n = which is identical to the free electron term of 
the n = 1 Landau level (n + l/2)h 2 s/2m. Nevertheless, 
the position is shifted by ~ 10 T between the two po- 
larizations. However, this can be understood when spin 
splitting is taken into account. The SWM prediction for 
the transitions are shown in Fig. Qb) together with the 
measured field splitting. It can be seen that there is in- 
deed no effect of electron-hole asymmetry for the 0—^1 
.PC-point transition (solid and broken lines). However, 
including the Zeeman term ig^i^B /2 with g = 2 the 
calculated splitting of the — > 1 transition, is compa- 
rable with the observed splitting. The Zeeman term is 
important here due to the very high magnetic field, and 
the fact that the — > 1 transition evolves very slowly 
with magnetic field so that a small energy splitting can 
generate a large field splitting in the data. 

Finally, we note that the Landau level energy spec- 
trum of graphene can be derived from the SWM Hamil- 
tonian simply by setting all the inter-layer coupling pa- 



rameters 71, ..,75 = 0. The analytic solution of this sim- 
plified Hamiltonian is nothing other than Eq.Q. This 
implies that the electron-hole asymmetry observed in 
the cyclotron resonance of exfoliated graphene [53] also 
originates from the neglected free electron kinetic energy 
terms. 
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